Supersymmetric Higgs pair discovery 
prospects at hadron colliders f 
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Abstract 

We perform a detailed study of the potential of hadron colliders in the search 
for the pair production of neutral Higgs bosons in the framework of the Minimal 
Supersymmetric Standard Model. The important role of squark loop contributions 
to the signal is emphasised. If the signal is sufficiently enhanced by these 
contributions, it could even be observable at the next run of the upgraded Tevatron 
collider in the near future. At the LHC the pair production of light and heavy Higgs 
bosons might be detectable simultaneously. 
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1. Introduction 

The search for Higgs bosons is one of the most 
important tasks for experiments at present and 
future high energy colliders ||l|. In particular, 
the Tevatron will soon start its next collider run, 
with slightly increased beam energy and greatly 
increased luminosity; a few years later experiments 
at the LHC will commence taking data. 

We study the production of two neutral Higgs 
bosons in gluon _fusion, followed by the decays of 
each boson into bb pairs. We focus on the final states 
where both Higgs bosons have (nearly) the same 
mass, since the resulting kinematical constraint 
helps to reduce the background. The SM cross 
section M is too small to be useful. However, the 
scalar sector of the SM suffers from well-known 
naturalness problems. These can be cured by 
introducing supersymmetry. Here we concentrate 
on the simplest potentially realistic supersymmetric 
model, the minimal supersymmetric standard 
model (MSSM) . Several effects can greatly enhance 
the Higgs pair production cross section in the 
MSSM with respect to the SM: 
1) If tan/3 3> 1, the Yukawa coupling of the b- 
quark is enhanced by a factor ~ tan/3 with respect 
to its SM value. It thus becomes comparable 



to the top quark Yukawa coupling for tan/3 ~ 
mt{mt) / m\,{mt) ~ 60, which is possible in most 
realizations of the MSSM. For Higgs boson masses 
around 100 GeV the squared 6-loop contribution 
then exceeds the i-loop contribution, which is 
suppressed by the large mass of the top quark, by 
a factor - 15 [|. 

2) For some region of parameter space {mA ~ 
300 GeV, tan/3 ^ 4) the branching ratio for H — > 
hh decays is sizable, h pair production through 
resonant H exchange is then enhanced by a factor 
{gMwIXtTH? -- 100 [|. 

3) Contributions from loops involving 6 or i squarks 
can exceed those from b and t quark loops by 
more than two orders of magnitude Q. This 
enhancement can occur for all values of niA and 
tan/3, but requires a fairly light squark mass 
eigenstate {ti or 6i), as well as large trilinear Higgs- 
squark-squark couplings. 

2. Monte Carlo Simulation 

In order to study the observability of the signal for 
Higgs pair production in the 46 final state, we have 
written MC generators for complete sets of signal 
as well as background processes. These generators 
were designed as new external user processes for the 
PYTHIA 5.7/JETSET 7.4 package §. 



We used the CompHEP package [|j to generate 
background events on the parton level. 

For both signal and background, the effects of 
initial and final state radiation, hadronization (in 
the string model) , as well as decay of the 6-flavoured 
hadrons have been taken into account. ( see |0]). 

3. Signal and Background Study 

We have calculated squark loop contributions to the 
pair production of two neutral Higgs bosons. If CP 
is conserved, squark loops contribute only if the two 
produced Higgs bosons have identical CP quantum 
numbers. We gave complete analytical expressions 
that allow the evaluation of these contributions 
(for details see |Q]). The Feynman diagrams 
contributing to the gg — > hh, HH, hH, and AA 
processes are presented in Fig. 1^. We take equal 
soft breaking contributions to diagonal entries of 
the stop and sbottom mass matrices {rrij: = m^ — 
m^ = rriq), as well as equal trilinear soft breaking 
parameters in the stop and sbottom sectors {At — 
Ab = Aq). We fix the running masses of the top 
and bottom quarks to mt{mt) — 165 GeV and 
mb{mb) = 4.2 GeV, respectively. This leaves us 
with a total of 5 free parameters, which determine 
our signal cross sections: ttia, tan/3, rriq, Aq and 
the supersymmetric higgsino mass parameter /i. 
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Figure 1. Feynman diagrams for hh, HH, hH, and 
AA Higgs boson pair production, -ffi(j) = h,H for 
i{j) = 1,2 respectively, qk(i) = qi,q2 for k{l) = 1,2. 



This parameter space is subject to experimental 
constraints 0, especially from the unsuccessful 
searches for Higgs bosons at LEP. We also demand 
that the masses of the lighter physical stop 
and sbottom exceed 90 GeV, which follows from 
squark searches at LEP. We also require that the 
contribution from stop and sbottom loops to the 
electroweak p-parameter satisfies 6p^^ < 0.0017. 
Finally, we only consider values of Aq and /i in the 
range \Aq\, |/i| < Sto^; this is necessary to avoid 
the breaking of electric charge and colour in the 
absolute minimum of the scalar potential. 



There are 6 different channels for producing 
two neutral Higgs bosons in the MSSM: HH, hh, 
AA, Hh, HA and hA. Often, several channels 
contribute to a given signal even after cuts have 
been applied, once the experimental resolution has 
been taken into account. The reason is that often 
two Higgs bosons are essentially degenerate in mass, 
especially for high tan/3. In our analysis we have 
combined contributions from different production 
channels assuming a Gaussian distribution for the 
reconstructed Higgs boson mass. We start with 
the diagonal process {hh, HH or AA production), 
which gives the best signal significance, and then 
add all other contributions to the "search window" . 
In order to give an idea of the signal rate for 
negligible squark loop contributions, we present in 
Fig. H contours of constant total signal cross section 
in fb in the (7nA,tan/3) plane. 
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Figure 2. Contours of constant cross section (in fb) for 
combined Higgs pair production channels, for the case 
of negligible squark loop contributions for the Tevatron 
(a) and the LHC (b). 

The total cross section is about 200 times higher 
at the LHC than at the Tevatron. Given an 
integrated luminosity of 100 fb~^, we expect well 
over 1,000 Higgs pair events at the LHC for all 
combinations of itia and tan/3. In contrast, if squark 
loop contributions are indeed small, at the Tevatron 
the raw signal rate is often too small to give a 
positive signal even at 25 fb^^ luminocity. 

In order to decide whether a Higgs pair cross 
section leads to a detectable signal, we have to 
compute the background rate. To suppress "fake" 
backgrounds, we require that all four 6-jets be 



tagged as such. The total cross sections for the two 
main irreducible backgrounds for the basic parton- 
level acceptance cuts pr > 25 GeV, ARjj > 0.5 
for the Tevatron (LHC) is 1.5 (59) pb for Zbb 
production and 2.6 (330) pb for bbbb production. 

The cross sections for the most important "fake" 
backgrounds for the Tevatron(LHC) are 3.1 (19.1] 
pb for Wbb [Q = A%) and 1.6 (164) nb for jjbb 
{Q = Mfji). Since the mis-tag probabihty of light 
quark and gluon jets is expected to be ^ 1% 
|g|,|l2l, after 5-tagging these "fake" backgrounds 
are much smaller than the irreducible backgrounds 
listed above and we therefore ignore them. 

One can see that irreducible backgrounds are 
clearly far larger than the signal. A more elaborate 
set of cuts is thus necessary. 

As already noted, we require all four 5-jets 
to be tagged. A realistic description of the b- 
tagging efficiency is therefore very important. In 
case of the Tevatron, we use the projected 5-tagging 
efficiency of the upgraded D0 detector g] and CMS 
collaboration iQ. We assume that 6-jets can be 
tagged only for pseudorapidity \ri},\ < 2 by both 
Tevatron and LHC experiments. 

We constructed the following kinematical 
variables and respective set of cuts for an efficient 
extraction of the signal: 

1) Reconstructed Higgs boson mass, Mh'- we 
chose the pairing that gives the smallest difference 
between the invariant masses of the two pairs: 
Mh = [Mb-^i,2+Mi,,ji,^]/2. After resolution smearing, 
the distribution in Mh for the signal can be 
described by a Gaussian with width a ~ ^/Mh 
(in GeV units). The search window is defined as: 
0.9m// in — 1.5cr < Mh < 0.9mH,in + 1.5(t. 

2) Mass difference between the masses of the two 
pairs (small for signal): iMbjf,^ — Mb^tJ < 2cr. 

3) The angles in the transverse plane between the 
two jets in each pair should be large, while the 
two transverse opening angles therefore tend to be 
correlated: 

A(/)6i,b2, A(t>b3M > 1, \A4>biM - A't'baMl < 1- 

4) All four 6-jets in the signal are fairly hard. We 
applied cuts on the softest and hardest of these jets, 
with transverse momenta pT,min and PT,max'- 
TEV ; pT,min > Mh/8 + 1.25a; pT,max > Mh/8 + 2a 
LHC : pT,min > Mh/4:; pT,-max > Mh/4: + 2a 

5) The 46 invariant mass: the signal distribution for 
this variable is concentrated around the invariant 
mass of the Higgs pair. This quantity has been 
shown to be useful for disentangling quark and 
squark loop contributions ji): M^h > 1-9Mh — 3cr. 
The efficiency of these cuts applied plus 46-tagging 
for several input (search) Higgs boson masses is 
listed in the following table for the Tevatron and 
LHC. 



The background efficiency refers to the cross 
section defined through the basic acceptance cuts 
(pt(6) > 25 GeV for all four 6 (anti-)quarks, and 
jet separation ARjj > 0.5 for all jet pairs). 
TEVATRON: 



mn.in [GeV 




120 


160 


200 


Esignal |%J 

ezbb [%] 
^bbbb [% 




2.10 
.187 
.137 


2.74 

.0935 

.0318 


3.30 

.0314 

.0072 


bbbb + Zbb 

Zbb 

bbbb 


# events 

for 

25 fb~i 


101 
11.1 
89.9 


26.3 

5.5 

20.8 


6.6 
1.9 

4.7 


signal [fbJ Br 
signal [fb] -Br 


95% c.l. 

5(7 


37.4 
76.9 


14.6 
30.0 


7.1 
19.4 


LHC: 


mH,in [GeV 




120 


160 


200 


^signal [%J 

ezbb [%] 

(■bbbb [% 


46 tag 


.34 

.0263 

.0142 


.90 

.0190 

.0112 


1.38 

.0081 

.0071 


bbbb + Zbb 
Zbb 

bbbb 


# events 

for 

100 fb-i 


4900 

240 

4660 


3863 

174 

3689 


2419 
73.6 
2345 


signal fb -Br 
signal fb -Br 


95% c.l. 
5cr 


570 
1426 


171 
427 


70.1 
175 



Table 1. Signal and background efficiencies and 

minimal cross sections for a 95% c.l. exclusion limit on, 
as well as a 5(7 discovery of, Higgs boson pair production 

This table also contains results for the minimal 
total signal cross section times branching ratio 
needed to exclude Higgs boson pair production at 
the 95% C.L, as well as the minimal total cross 
section times branching ratio required to claim a 
5(7 discovery of Higgs boson pair production in the 
46 final state. We give these critical cross sections 
for two values of the integrated luminosity at the 
Tevatron, characteristic for the upcoming Run II 
and for the final luminosity at the end of the TeV33 
run, respectively. 

Systematic uncertainties are a concern, espe- 
cially at the LHC, where the large signal rate can 
lead to a very small signal-to-background ratio if 
the significance is defined using statistical errors 
only. We assign a systematic uncertainty of 2% on 
the background estimate, as obtained by extrapola- 
tion from the side bins. We thus require a minimal 
signal-to-background ratio of 0.04 for the 95% c.l. 
exclusion limit, and 0.1 for the 5(7 discovery cross 
section. This requirement in fact fixes the critical 
cross sections at the LHC for iriH.in < 180 GeV. 

4. Potential of Hadron Colliders for Higgs 
Pair Search 

By comparing the results of Table IT] and Fig. pla, it 
becomes clear that in the absence of sizable squark 
loop contributions to the signal cross section, the 
potential of Tevatron experiments for this search 



is essentially nil. In contrast, some parts of the 
(rriy^i, tan/3) plane can be covered at the LHC even 
if squark loop contributions are negligible. For 
this pessimistic assumption of negligible squark loop 
contributions, LHC experiments might discover a 
5a signal if tan/3 is large ( ^ 50), and can at least 
exclude some regions of parameter space where tan/3 
is small ( ^ 2.5). 

In order to illustrate the possible importance 
of squark loop contributions, we performed various 
Monte Carlo searches of the three-dimensional 
parameter space (m^, Aq, fi). We believe that 
our procedure should reproduce the maximal cross 
section to within a factor of 2 or so. 




Figure 3. 95% c.l. exclusion and 5a discovery regions 
for Higgs pair production at the Tevatron (25 fb~^) (a) 
and LHC (100 fb~^) (b), for "maximized" squark loop 
contributions. The light grey contour in (a) shows the 
region where a, > 5a signal should be detectable at the 
Tevatron with just 2 fb~^ of data. 

The results are presented in Fig. 0, which shows 
the regions that can be probed with 2 and 25 fb^^ 
of data at the Tevatron (a), and with 100 fb~^ of 
data at the LHC (b). We see that now virtually 
the entire part of the (rn^, tan/3) plane will give 
a > 5(7 signal at the LHC. Moreover, the entire 
region niA < 200 GeV, and most of the region with 
niA < 300 GeV, can be probed at the Tevatron 
with 25 fb~^ of data. Perhaps the most surprising, 
and encouraging, result is that a substantial region 
of parameter space will give a > 5cr signal at the 



Tevatron already with 2 fb^^ of data! This is the 
first time that such a robust signal for Higgs boson 
production at the next run of the Tevatron collider 
has been suggested. 

5. Summary and Conclusions 

The main outcome of this analysis consists in values 
of the minimal total signal cross section times 
branching ratio required for a 5a observation of 
the signal, as well as for placing 95% c.l. exclusion 
limits, at both the Tevatron and the LHC. 

In the absence of substantial squark loop con- 
tributions, the prospects for Tevatron experiments 
appear to be dim. LHC experiments can then only 
probe scenarios with niA ~ 300 GeV and either very 
large or quite small values of tan/3. 

On the other hand, if squark loop contributions 
are nearly maximal, and if it is possible to construct 
an efficient trigger for events containing 4 6-jets with 
(pt) ^ 50 GeV, LHC experiments should find a 
signal for hh production for practically all allowed 
combinations of rriA and tan/3; HH production 
(augmented by nearly degenerate modes) should 
be visible for most scenarios with ttih < 
2mt. Moreover, with 25 fb~^ of data, Tevatron 
experiments would be sensitive to most of the region 
with TTiA < 300 GeV; if tan/3 is large, even scenarios 
with uiA > 500 GeV might be detectable. Our 
most exciting result is that a significant region of 
parameter space with niA ^ 250 GeV should be 
accessible already at the next run of the Tevatron 
collider, which is projected to collect 2 fb~^ of data. 
This seems to be the most robust signal for the 
production of MSSM Higgs bosons at the Tevatron 
that has been suggested so far. 
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